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1. Introduction 

With twisted mass lattice QCD(tmQCD) [JIJ], it became possible to perform lattice simulations 
with Wilson type quarks at small quark masses. This is because of the fact that the Wilson Dirac 
operator with a twisted mass term is free of unphysical zero eigenmodes as long as the quark mass 
and the twist angle are nonzero. Another extremely useful property of tmQCD is that when the 
twist angle equals 5 (maximal twist), discretization errors of 0{a 2k+l );k = 0,1,2,... of correla- 
tion functions of gauge invariant, multiplicatively renormalizable, multi-local operators are either 
automatically absent or could be removed by taking a simple averaging procedure^]. This is par- 
ticularly simple in the case of extracting masses from 2-point correlators, since the improvement 
will be automatic and no averaging is needed. These features make tmQCD at maximal twist an 
attractive method to study the hadronic spectrum. An important element in this program is the 
definition of maximal twist and how the uncertainty in this definition could affect the expected 
reduction of discretization errors. Initial quenched simulations [||] used the condition of maximal 
twist from the vanishing of the pion mass as computed in the untwisted action (Wilson maximal 
twist). Improved results were obtained using Wilson maximal twist at not so light quark 

masses. However, for very light quark masses (close to the interesting physical up and down 
quark masses), undesired discretization effects became evident[Q][Q]. In addition tmQCD chiral 
perturbation theory [jv[] [J8|] as well as a study of discretization errors in tmQCD a la Symanzik [§] 
argued for the need for an optimal determination of the maximal twist condition in the light quark 



mass regime. Improved definitions of maximal twist Ql lfl [ |12| ] have been used in recent simulations 



J^][|l(]][|i~ri], using PCAC or the vanishing of wrong parity matrix elements(parity maximal twist). 



For more discussion, see Ql3p. 

In this work, results for the quenched hadron spectrum from a degenerate doublet of up and 
down quarks are presented for various quark masses and at different lattice spacings. Simulations 
were performed implementing both the Wilson maximal twist and parity maximal twist definitions 
and results are compared. 

2. Action at Maximal Twist 

The action is given by 

S[yr, w,U,p]=S F [y, % U] +S g [f5,U], (2. 1) 

where S g is the standard Wilson plaquette action for the gauge field at coupling j3, while Sf is the 
tmQCD action for a degenerate doublet y of the up and down quarks. In the twisted basis 

S F [yr, % U] = a 4 £ \ff(x) (M + ipft t 3 + %) y(x) , (2.2) 

X 

where 

^w = ^I(7vV v + 7vV;-«V;v v ), 

Z V 

V v Y(x) = U v (x)y(x + av) - \jf(x), V* v \j/(x) = Xjf(x) - Ul(x - av)\ff(x - av) . 
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#sites 


Configurations 


aMcritical 


a\l q 


physical 

am' q J 


5.85 


20 3 x 40 


1000 


-0.9071 


0.0376 
0.0188 
0.012527 
0.00627 


~ m s 
~ m s /2 
~ m. s ./3 
~ m . s /6 






600 


-0.8965 
-0.0971 
-0.9110 
-0.9150 


0.0376 
0.0188 
0.01252 
0.00627 


~ m s 
~ m s /2 
~ m. s /3 
~ m s /6 


6.0 


20 3 x 48 


1000 


-0.8135 


0.030 
0.015 
0.010 
0.005 


~ m s 
~ m s /2 
~ ra. s /3 
~ m s /6 






600 


-0.8110 
-0.8170 
-0.8195 
-0.8210 


0.030 
0.015 
0.010 
0.005 


~ m . v 
~ m^/2 
~ m. v /3 
~ m^/6 


6.2 


28 3 x 56 


200 


-0.7363 


0.02165 
0.01083 
0.00722 
0.00361 


~ m . v 
~ m s /2 
~ m s /3 
~ mi/6 



Table 1: Parameters of the simulations. The bare quark mass apt q values were chosen such that the physical 
quark mass am,q hys,cal corresponds to specific fractions of the strange quark mass. For Wilson maximal twist, 
the value of aM cr itical i s independent of the bare quark mass value. 



When the mass parameter M = M cr ni ca i, we are at maximal twist and the bare quark mass equals 
ajl q . Simulations are performed with M critica i defined by the value that corresponds to a massless 
pion in the standard Wilson action (Wilson maximal twist) [||], and with M cr jti C al defined by the 
value when the wrong parity mixing between the physical vector and pseudoscalar vanishes (parity 
maximal twist), i.e. 

£<V- (x,t)P + (0,0) > =0. (2.3) 

X 



3. Simulations 

In Table [I], the simulation parameters are listed. The values of M cr iti ca i for the parity definition 



of maximal twist were tuned at each \i q such that the condition in Eq.jlJjis satisfied. Standard local 



interpolating fields for the charged pion, the charged p , J p = (^) ± baryons and J p = (1)^ baryons 
are used to obtain the ground state masses[||]. Using the PCAC relation, the charged pion decay 
constant can be obtained from 

= % <0|P»|. (3.1) 
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P=5.85 



P=6.()0 



P=6.2 

i r 




Figure 1: The charged pion mass squared (in lattice units) as a function of the bare quark mass a[i q . 



4. Results 

In Figure [I], the charged pion mass squared is plotted as a function of the bare quark mass. The 
parity definition of maximal twist shows the expected linear dependence at light quark masses with 
am n approximately vanishing in the chiral limit. A better scaling behavior of f % is also observed 
with the parity definition of maximal twist as shown in Figure ||, while the scaling of the rho mass 
is equally good in both definitions if one takes into account the large statistical errors. For the 
rho mass, results extracted using the 4 th component of the tensor operator are also shown which 
were found, in some cases, to have small statistical errors. In Figure ||, results for the Nucleon and 
Delta masses are shown for both definitions of maximal twist. These are about 10% higher than 
the physical values, a result that might be expected due to quenching. 

The twisted mass action includes flavor breaking effects due to the Wilson term. These could 
be seen as a mass splitting among isospin multiplets. In case of the A multiplet, there are no 
disconnected diagrams and one can study the mass splitting among the four A states. On the 
average, for a degenerate doublet of the up and down quarks, the up quark propagator U (x,y) and 
the down quark propagator D(x, y) are not equal but are related by [^] 



U{x,y) = 7 5 £>(y,x) + 75- 



(4.1) 



Using this relation, one expects the following mass relation among the A states, 



(M A ++ = M A - =M A ++,-) ^ (M A + =M A o =M A+ ,o). 



(4.2) 



In Figure f|, the measured mass splittings are shown for /3 = 6.0. From these results, one concludes 
that there could be a hint of flavor breaking, however, these are hidden by large statistical errors. 
The flavor breaking effects between the charged and neutral pion need the computation of discon- 



nected diagrams which has been calculated in [14] [15]. A recent study of tmQCD with a strange 



quark also includes some discussion of flavor breaking] 16]. 
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Figure 2: Scaling of the charged pion decay constant and the charged rho mass. 
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Figure 3: Baryon mass for (^) + and (|) + states as a function of the pion mass squared. The physical 
nucleon N(939) and A(1232) masses correspond to the lattice spacings in the legends. 



5. Conclusions 

Improved results for the hadron spectrum could be obtained from tmQCD and light quark 
masses. The parity definition of maximal twist was found to have a better scaling and chiral be- 
haviour than the Wilson definition of maximal twist. 
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Figure 4: A mass splittings. 
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